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ABSTRACT 


As  ownship  moves  through  a  sea  of  targets  moving  at  eonstant  veloeity  in  random 
direetions,  it  will  oeeasionally  meet  some  of  them  when  the  intervening  distanee  beeomes 
suffieiently  small.  The  very  faet  of  meeting  is  signifieant  in  estimating  the  properties  of  a 
met  target.  We  derive  joint  distributions  for  the  properties  of  met  targets,  and  show  how 
to  sample  those  properties  in  Monte  Carlo  simulations.  We  also  introduce  “instant”  target 
motion  analysis  in  the  form  of  various  conditional  distributions  of  target  properties,  given 
observations  made  at  the  time  of  meeting. 
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1. 


INTRODUCTION 


Equation  Section  IWe  consider  situations  where  a  certain  vehicle,  called 
“ownship”,  moves  at  constant  velocity  through  an  unbounded,  two-dimensional  medium 
full  of  other  vehicles  that  are  themselves  moving  at  constant  velocities  in  independent 
random  directions.  It  sometimes  happens  that  ownship  meets  one  of  these  other  vehicles 
(a  “target”),  by  which  we  mean  that  the  distance  between  them  decreases  to  some  critical 
value.  The  details  of  such  meetings  can  be  important  to  both  parties — a  collision  may  be 
imminent,  or  the  target  may  be  about  to  commit  a  hostile  act  that  requires  immediate 
action  by  ownship,  with  the  action  depending  on  the  details.  The  detection  mechanism 
strongly  influences  what  can  be  said  about  those  details.  If  detection  is  due  to  an  active 
sensor  like  radar,  then  an  accurate  target  location  may  be  immediately  available,  but  there 
are  also  passive  sensors  that  can  sometimes  announce  the  fact  of  detection  without 
providing  much  detail,  especially  detail  about  the  distance  to  the  target. 

This  subject  is  particularly  important  for  ships,  and  so,  even  though  there  is 
nothing  particularly  nautical  about  the  rest  of  this  paper,  we  will  refer  to  the  vehicles  as 
ships  from  here  on.  We  will  also  use  nautical  terms  such  as  “range”  (distance  to  the 
target),  “course”,  and  “bearing”,  but  all  angles  are  measured  counterclockwise  in  radians 
from  East,  the  conventional  mathematical  system. 

Even  though  all  ships  are  assumed  to  be  moving  uniformly  at  random  a  priori, 
that  statement  is  not  true  of  targets  that  are  met.  It  is  often  far  from  being  true,  but 
nonetheless  it  will  be  tempting  for  an  analyst  to  suppose  that  the  course  or  bearing  of  a 
met  target  is  (say)  uniformly  distributed  over  some  plausible  range.  In  this  paper,  we  will 
derive  nonuniform  distributions  for  such  quantities  that  follow  directly  from  the  a  priori 
assumptions  and  Bayes’  theorem.  Monte  Carlo  simulations  (e.g.,  Eorman,  2005) 
sometimes  have  ownship  meeting  a  succession  of  targets  as  it  moves,  with  each  met 
target  being  assigned  properties  such  as  range,  bearing,  speed,  and  course.  The 
aforementioned  distributions  can  be  used  to  set  these  properties  realistically.  Section  6.3 
includes  some  sampling  functions. 

The  careful  distinction  between  targets  in  general  and  targets  that  have  been  met 
goes  back  to  at  least  World  War  Two.  Let  b  be  the  (relative)  bearing  from  ownship  to  a 
met  target  at  the  moment  of  meeting,  by  which  we  mean  the  angle  from  ownship ’s 
velocity  vector  to  the  line-of-sight  to  the  target  (Eigure  1).  Meetings  are  more  likely  to 
occur  with  small  bearings  (bow  sightings)  than  large  ones  (stem  sightings);  indeed,  stem 
sightings  are  impossible  if  ownship  has  a  higher  speed  than  the  target.  It  follows  that 
moving  sensors  desiring  to  detect  targets  as  early  as  possible  should  spend  more  of  their 
time  looking  near  the  bow  than  in  other  directions.  The  density  function  of  b,  given 
detection  of  a  target  with  known  speed,  was  derived  in  World  War  Two  (Koopman, 
1956)  in  the  process  of  dealing  quantitatively  with  the  issue  of  where  lookouts  should 
concentrate  their  effort.  This  paper  generalizes  that  work. 

Although  we  anticipate  that  the  main  application  of  our  results  will  be  in  setting 
the  parameters  of  newly  met  targets,  there  are  also  potential  applications  to  Bayesian 
inference.  The  most  notable  of  these  is  the  determination  of  the  target’s  range  from 
bearing  measurements,  sometimes  known  as  target  motion  analysis  (TMA).  It  is  well 
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known  (Wagner,  Mylander,  &  Sanders,  1999,  Theorem  1 1.4)  that  range  eannot  be  exaetly 
inferred  from  bearing  measurements,  no  matter  how  numerous  and  accurate,  unless 
ownship  maneuvers  (that  is,  unless  ownship  departs  from  the  constant  course  and  speed 
assumptions  that  we  make  throughout  this  paper).  However,  that  fact  does  not  prevent  the 
application  of  Bayes  theorem — an  activity  that  we  dub  “instant”  TMA  to  distinguish  it 
from  single-leg  TMA  where  ownship  does  not  turn,  but  makes  a  sequence  of 
measurements  separated  in  time  (e.g.,  Ganesh  &  Baylog,  2004). 

If  ownship  detects  a  target  using  passive  sensors,  it  is  likely  that  b  and  perhaps  its 
time  derivative  z  =  dbldt  (the  bearing  rate)  will  be  known.  Those  known  quantities  can 
then  be  used  in  instant  TMA.  We  will  not  consider  errors  in  measuring  either  of  those 
quantities.  Doing  so  would  be  essential  if  we  would  otherwise  obtain  exact  estimates  of 
operational  quantities,  but  that  is  not  the  case  when  one  is  merely  applying  Bayes’ 
theorem  to  data  measured  or  discovered  at  one  point  in  time.  Intuition  is  that  the  target 
must  be  nearby  if  the  bearing  rate  is  high.  This  turns  out  to  be  correct,  but,  as  will  be  seen 
in  Section  7,  considerable  uncertainty  about  the  target’s  range  will  remain  even  when  b 
and  z  are  known  exactly.  Instant  TMA  typically  provides  only  a  slight  refinement  in 
knowledge  of  operational  quantities  such  as  the  target’s  range. 

Problems  where  the  speed  and  detection  range  of  the  met  target  are  either  known 
or  irrelevant  are  taken  up  in  Section  4.  Such  problems  include  Koopman’s  derivation  of 
the  bearing  distribution  referred  to  above.  For  other  questions,  an  application  of  Bayes’ 
theorem  will  require  a  joint  a  priori  distribution  for  the  target’s  speed  and  detection 
range.  A  useful  class  of  such  distributions,  based  on  the  passive  radar/sonar  equation,  is 
introduced  in  Section  5. 

The  reader  may  wish  to  have  available  the  Excel™  workbook  Meetings.xls,  which 
contains  Visual  Basic  code  for  some  of  the  aforementioned  distribution  functions,  and 
also  illustrates  their  application.  That  workbook  can  be  found  among  the  downloads  at 
http  ://facultv.nps  .edu/  awashbum/. 
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2. 


BASIC  DEFINITIONS  AND  RELATIONSHIPS 


Figure  1  shows  two  ships,  each  moving  in  a  straight  line  with  constant  velocity. 
Without  loss  of  generality,  we  assume  ownship  is  proceeding  to  the  right.  The  target  is  at 
bearing  b  and  range  r,  proceeding  at  speed  u.  The  target’s  course  is  c,  relative  to 
ownship ’s  course,  or  \\i,  relative  to  the  line-of-sight  (LOS)  from  ownship  to  the  target 
(V  =  c-b). 

Figure  2  shows  the  same  situation  in  a  coordinate  system  where  ownship  remains 
stationary.  The  velocity  of  the  target  relative  to  ownship  is  the  vector  difference  between 
target  velocity  and  ownship  velocity,  shown  in  both  polar  and  Cartesian  forms.  The 
target’s  relative  speed  is  w  and  its  relative  course  is  (j),  or  alternatively,  the  target’s 
velocity  can  be  resolved  into  two  components:  one  (x)  within  the  LOS  and  directed 
toward  ownship,  and  the  other  (y)  normal  to  the  LOS  and  directed  counterclockwise.  The 
formulas  relevant  to  our  analysis  are: 

x  =  v  cos(h)  -  u  cos(c  -  b) 

y  =  vsm{b)  +  usm{c-b)  (1.1) 

w  =  Ju^  +v^  -2mvcos(c) 


Figure  1:  Ownship  and  target  moving  at  constant  velocity  in  two  dimensions. 

The  range,  bearing,  speed,  and  course  of  the  target  are  all  assumed  to  be  random 
variables,  and  uppercase  versions  of  r,  b,  u,  and  c  will  be  used  to  represent  them.  The 
general  idea  is  that  the  two-dimensional  medium  is  initially  assumed  to  be 
homogeneously  populated  with  targets  located  according  to  some  constant  Poisson 
density.  Each  of  these  targets  is  then  assigned  parameters  {R,  U,  C),  independently  of  the 
others.  Specifically,  R  and  U  are  assigned  according  to  some  density  g{r,u)  that  will  be 
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left  indefinite  for  the  moment,  and  C  is  independently  assigned  to  be  uniform  over  a  full 
eirele.  Eaeh  target  moves  aeeording  to  its  assigned  speed  and  eourse,  and  ownship  moves 
to  the  right  at  speed  v.  Meetings  oeeasionally  happen  when  the  distanee  between  ownship 
and  some  target  decreases  to  the  detection  range  assigned  to  that  target  (we  ignore  any 
initial  detections  that  happen  when  R  is  larger  than  the  initial  distance  to  the  target).  Each 
of  these  meetings  is  characterized  by  parameters  {R,  B,  U,  and  C),  and  we  wish  to 
describe  various  distributions  and  conditional  distributions  involving  those  four  random 
variables.  The  joint  distribution  of  the  four  is  the  most  fundamental,  and  the  subject  of  the 
next  section. 


Figure  2:  Illustrating  target  velocity  relative  to  ownship,  resolved  into  components  x 
directly  toward  ownship  and  j  across  the  LOS,  Meetings  can  happen  only  when  x 

is  positive. 

The  assumption  that  C  is  a  priori  uniform  over  a  full  circle  is  significant,  and 
could  be  false  if,  for  example,  the  targets  were  all  merchant  ships  moving  within  a  traffic 
lane.  The  uniformity  assumption  is  probably  more  realistic  for  military  ships  than  for 
merchant  ships. 

Eor  brevity,  we  will,  in  most  cases,  use  the  notation XEh.w.c)  for  the  joint  density, 
in  place  of  the  more  precise  noidXion  fR_B,u,dr’b,u,c),  trusting  the  name  of  the  dummy 
variable  to  cue  the  name  of  the  random  variable  under  discussion,  and  likewise  for 
discovery  rates  such  as  X{r,b,u,c).  Similarly,  y(-^[y)  is  the  density  of  random  variable  X 
whenever  the  random  variable  Y is  given  to  bey. 
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3. 


THE  JOINT  DENSITY  OF  (R,  B,  U,  C)  UPON  MEETING 


First  consider  some  elass  of  targets,  all  of  which  have  speed  u  and  detection  range 
r,  and  assume  that  these  targets  are  initially  loeated  in  a  two-dimensional  Poisson  field 
relative  to  ownship,  with  density  h  targets  per  unit  area.  All  targets  initially  eloser  than  r 
to  the  origin  (ownship’s  location)  are  eliminated  from  eonsideration,  since  they  are 
detected  immediately.  Eaeh  of  the  remaining  targets  is  independently  assigned  a  constant 
course  C  that  is  a  uniform  random  variable  in  the  interval  [-tt,  tt].  Should  one  of  these 
targets  be  subsequently  met,  it  will  also  have  a  bearing  B,  as  well  as  a  range  r  from 
ownship.  Following  Koopman’s  1956  argument,  the  rate  at  whieh  ownship  meets  sueh 
targets  with  speeifie  eourse  c  and  bearing  b  is 

A(r,b,u,c)  =  2rhx^ (db)(dc  /  Itt);  b,c  e  ,  (1-2) 

where  x  is  the  approaeh  veloeity  given  by  (1.1)  and  the  +  sign  means  “positive  part”;  that 
is,  is  X  if  X  >  0  or  0  if  X  <  0.  It  should  be  intuitively  elear  that  meetings  are  impossible  if 
X  <  0,  since  sueh  targets  are  moving  away  from  ownship,  rather  than  toward  it.  The 
produet  of  2r  with  x^  is  an  area  eoverage  rate. 

By  integrating  (1.2)  over  all  values  of  b  for  whieh  x  is  positive,  we  find  that 

A.{r,u,c)  =  2rhw{dc  / Itt)  =  2rh^jM^  +v^  -2uvcos(c)  {dcl2ny,  c  e  (1-3) 

where  w  is  the  relative  speed  from  (1.1).  If  we  now  integrate  over  all  values  of  c,  we 
obtain  the  rate  of  finding  targets  of  this  elass: 

Itt 

X{r,u)  =  2rh — [  +v^  -2Mveos(c)  dc  =  2rhS{u,v)  .  (1-4) 

0 

The  equivalent  speed  function  S{u,v)  is  a  complete  elliptie  integral  of  the  seeond 
kind.  It  represents  the  expected  relative  speed  between  two  targets,  with  speeds  u  and  v 
moving  in  uniformly  random  direetions.  The  function  X{r,u)  is  the  rate  parameter  for  a 
Poisson  proeess  of  target  diseovery.  If  w  =  0,  S{u,v)  reduees  to  v,  and  (1.4)  reduees  to  the 
familiar  idea  that  the  rate  of  diseovering  targets  is  the  rate  of  examining  area  (2rv) 
multiplied  by  the  density  of  targets  per  unit  area  (h). 

Now  suppose  that  eaeh  target  is  initially  assigned  a  detection  range  and  speed 
{R,U)  aeeording  to  the  density  funetion  g{r,u),  in  whieh  ease  the  diseovery  rate  for  targets 
with  speeifie  {dr,db,du,dc)  is  Z{r,b,u,c)g{r,u){dr){du)  .  By  integrating  over  b  and  c  as 
before,  the  diseovery  rate  for  targets  with  speeifie  range  and  speed  {dr,  du)  is 
X{r ,u)g{r ,u){dr){du) .  Let 

00 

Rg{u)  =  \rg{r,u)dr .  (1.5) 

0 
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Then  the  total  rate  of  diseovering  targets  is 


jiMr  ,u)g(r,u)(dr)(du)  =  2A| R^(u)S(u,v)(du)  =  2hA  .  (1-6) 

0  0  0 

The  quantity  2A  has  dimensions  of  area  per  unit  time,  and  ean  be  thought  of  as  an  overall 
rate  of  examining  area  for  targets. 

The  joint  density  funetion  of  {R,B,U,C)  for  a  diseovered  target  is  now  given  by  a 
ratio  of  rates; 


fir  b  u  c)=  ^  rgir,u)iv  eos(h)  -  u  eos(c  -  b)y 

,u,c)  ^  .  t  •  ; 

Note  that  h  eaneels  when  taking  this  ratio,  so  the  spatial  density  of  targets  is  immaterial. 
The  joint  density  (1.7)  ean  be  used  to  derive  various  other  quantities.  In  partieular,  the 
joint  density  of  iR,U)  for  a  met  target  is 

f{r,u)  =  X{r ,u)g{r ,u) l{2hA)  =  rg{r,u)S{u,v)IA  .  (1-8) 


Note  that y(r,M)  is  not  gir,u).  Due  to  the  faetor  r  in  (1.8),  met  targets  tend  to  have 
longer  deteetion  ranges  than  do  targets  in  general.  Fast  targets  also  tend  to  be  over¬ 
represented  beeause  Siu,  v)  inereases  with  u.  This  is  in  eontrast  to  the  observed  faet  that 
fast  fish  are  underrepresented  in  towed  sampling  nets,  presumably  beeause  fish  use  their 
speed  to  avoid  the  net  when  they  see  it  (Barkley,  1964),  rather  than  eontinuing  to  move  at 
eonstant  veloeity  as  we  assume  here. 

The  density  of  U upon  meeting  ean  be  obtaining  by  integrating  (1.8)  over  all  r: 


R  (u) Siu, v) 
fiu)  =  ^-^ - 


(1.9) 


In  faet,  the  joint  density  of  any  subset  of  ir,b,u,c)  ean  be  obtained  by  integrating  the  joint 
density  X^.h.M.c)  over  the  variables  not  in  the  subset.  Equation  (1.7)  is  the  joint  density, 
and  all  others  ean  be  derived  from  it. 

If  there  are  aetually  several  elasses  of  targets  indexed  by  i,  with  elass  i  having 
spatial  density  hi  and  initial  density  giir,u)  for  iR,U),  then  the  rate  of  diseovering  targets 
of  elass  i  is  2/i  A,  where 

00  00 

Aj.  =^Siu,v)i^rg^ir,u)du)dr ,  (1-10) 

0  0 

and  the  total  rate  of  diseovering  targets  is  D  =  ^  2/i.A .  .  If  the  diseovery  proeess  for  eaeh 

i 

elass  of  targets  is  an  independent  Poisson  proeess  with  rate  2/i,A/,  then  the  diseovery  of 
all  targets  is  a  Poisson  proeess  with  overall  rate  D.  Equivalently,  one  ean  imagine  an 
overall  Poisson  diseovery  proeess  with  rate  D  where  the  target  elass  is  determined  upon 
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discovery  to  be  i  with  probability  2A/A,/Z),  with  the  elass  determining  the  joint 
distribution  of  {R,B,U,C).  Here  we  have  used  results  about  the  classification  and 
superposition  (summing)  of  independent  Poisson  processes  (Heyman  &  Sobel,  1990; 
Ross,  2000).  We  note  these  multielass  possibilities  only  in  passing,  and  will  eonsider  only 
a  single  class  from  here  on.  Section  4  considers  quantities  that  do  not  involve  the  as  yet 
unspecified  prior  density  g(r,u).  Section  5  deals  with  inferences  that  require  specific 
assumptions  about  g(r,M). 
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4. 


DISTRIBUTIONS  INDEPENDENT  OF  g{r,u) 


4,1  JOINT  DENSITY  OF  (S,C),  GIVEN  U:  7(6, c|a) 

In  random  meetings,  the  target’s  bearing  B  and  eourse  C,  relative  to  the 
line-of-sight,  are  neither  independent  nor  uniform.  From  (1.2)  and  (1.4), 


f{b,c  \r,u)  =  f{r,b,u,c)/ fir, u)  =  - - 

4;ro  (m,v) 


(v  cos(6)  -  u  cos(c  -  b)y 
47rSiu,v) 


(1.11) 


Sinee  the  right-hand  side  of  (1.1 1)  is  independent  of  r,  the  eonditional  density  is  the  same 
thing  asfib,c\u).  Equivalently,  sinee  \|/  =  c-b. 


fib,y/\u)  = 


47iSiu,v) 


(v  cos(b)  -  u  cos((//)y 
47rSiu,v) 


(1.12) 


4,2  DENSITY  OF  B,  GIVEN  U:  fib\u) 

To  find  7(6 1 m),  we  need  to  integrate  the  density  (1.12)  over  all  \|/  for  whieh  is 
positive.  If  veos(6)  <  -u,  thenx^  =  0  for  all  \|/,  and  oonsequently7(6|M)  =  0.  If  voos(6)  >  u, 
then  =x  for  all  \|/,  and,  sinee  the  integral  of  oos(\|/)  over  a  full  eirele  is  0, 
fb\u)  =  veos(6)/(25'(M,v)).  There  remains  the  ease  where  |veos(0)|  <  u,  in  whieh  ease  there 
is  an  angle  \|/o  sueh  that  Meos(\|/o)  =  veos  (6)  and  0  <  \|/o  <  ir.  Sinee  >  0  if  and  only  if 
\|/o  <  x]/  <  271  -  \|/o,  by  symmetry 


K 

ivcosib)-ucosii//))cb// 


fib\u)  =  ^ 


47rSiu,v) 


(v  eos(6)(;r  -y/y  +  u  sin(^o )) 
2nSiu,v) 


Sinee  tt  -  \|/o  =  AreCos(-veos(6)/M),  and  sinee  Msin(\j/o)  = ,  the  above  can 
be  summarized  as 


fib\u)  = - 

2nSiu,v) 


0,  if  vcos(6)  <  -u 

<  ;rv  cos(6),  if  V  cos(6)  >  M  (1-13) 

V cos{b)ArcCos{-v cos(6) / w)  +  sju^  -v^  cos^(6),  if  \vcosib)  |<  u. 


This  is  the  density  first  determined  by  Koopman  (1956)  in  World  War  Two.  Targets 
cannot  appear  outside  the  “limiting  lines  of  approach”  (the  first  condition),  and  are  more 
likely  to  appear  over  the  bow  (near  6  =  0)  than  near  any  other  bearing.  Only  the  ratio  of 
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the  two  speeds  is  relevant.  Figure  3  shows  the  distribution  when  m/v  =  0.5,  1, 
and  2. 


Sinee  r  is  missing  from  the  right-hand  side  of  (1.11),  formula  (1.13)  is  also 
J{b\r,u);  the  additional  knowledge  of  the  detection  range  is  of  no  value  as  long  as  the 
target’s  speed  is  already  known. 


- u/v=0.5 

. u/v=1.0 

- u/v=2.0 


bearing  (radians) 


Figure  3:  Density  of  the  bearing  at  meeting  (B)  for  three  speed  ratios. 

4.3  DENSITY  OF  \|/,  GIVEN  B  AND  U:  fiMb,u) 

Given  that  a  target  has  been  detected  on  bearing  b,  this  is  the  density  of  the 
target’s  course  relative  to  the  line  of  sight,  which  is  of  interest  in  forecasting  the  target’s 
motion.  This  density  is  simply  the  ratio  of  (1.12)  to  (1.13).  Note  that  S{u,v)  cancels  in 
taking  the  ratio,  so  this  is  the  easiest  to  evaluate  of  any  of  the  densities  in  this  section.  If 
the  ratio  is  0/0,  the  interpretation  is  that  b  is  impossible. 

4.4  JOINT  DENSITY  OF  B  AND  T,  GIVEN  U: 

Our  interest  in  Y,  the  target’s  speed  across  the  LOS,  stems  from  the  random 
variable’s  close  relationship  to  Z,  the  bearing  rate,  specifically  Y=RZ.  Since  Z  is 
passively  observable,  we  can  potentially  condition  on  it,  as  well  as  on  B,  when  applying 
Bayes’  theorem. 

Since  Tis  a  function  of  both  B  and  \|/  according  to  (1 . 1),  this  section  is  an  exercise 
in  calculating  the  density  of  a  scalar  function  of  two  random  variables,  complicated  by 
the  fact  that  the  function  involved  is  not  mono  tonic.  In  general,  if  the  angles  for  which 
y(h,\i/)  =  y  are  \|/i,  \i/2,  . . . ,  then/(h, y\u)  =  E,  f{b,i//.\u)/\  Y'{y/^)  I ,  where  ^(xj/i)  denotes 
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the  derivative  of  Y  with  respeet  to  \j/  at  the  f'  angle.  In  the  present  application,  the  sum 
may  include  0,  1,  or  2  terms.  From  (1.1),  the  required  derivatives  are  all  F(\|/i)  =  MCOs(\|/i). 

There  are  multiple  cases.  Refer  to  Figure  4.  If  vcos(b)  <  -u,  i\iQnj{b,y\u)  is  0  for 
all  y  because  b  is  impossible  (x  is  negative  regardless  of  the  target’s  course).  Otherwise, 
let  5=y-vsin(h),  so  that  the  equation  Y{b,\y)=y  is  equivalent  to  the  equation 
5  =  Msin(\|/).  If  |5|  >  u,  then  j{b,y\u)  =  0  because  there  are  no  solutions  to  the  equation. 
Otherwise  there  are  two  solutions,  namely  \|/i  =  ArcSin(5/M)  and  \|/2  =  7t:-\|/i.  For  these 

two  solutions  we  have  MCOs(\j/i)=  -d^  and  mcos(\|/2)  =  -mcos(\|/i),  with  each  solution 
qualifying  for  inclusion  in  the  sum  if  and  only  if  the  closure  rate  x  at  that  solution  is 
nonnegative.  If  vcos(h)  >  mcos(\)/i),  then  both  qualify,  and  the  sum  of  the  two  joint 
densities  is  j{b,\yi\u)+j{b,\y2\u)  =  2vcos{b)l{AiiS{u,v)).  If  mcos(\|/2)  <  vcos(h)  <  mcos(\|/i), 

then  only  \|/2  qualifies,  with  joint  density  / {b, I  w)  =  (v cos(h)  +  -sju^  -S^  )l{A7rS{u,v))  . 
If  vcos(h)  <  mcos(\|/2),  then  neither  solution  qualifies.  In  all  cases,  the  joint  density  must 
be  divided  by  |Y'(\|/i)|,  which  is  ,  before  being  summed. 


y 


i/sin(ib) 


Figure  4:  Velocity  across  the  line  of  sight  (y)  on  the  vertical  axis,  and  closing  velocity 
along  the  line  of  sight  (x)  on  the  horizontal  axis.  The  angles  \|/i  and  \|/2  both  qualify  in 
this  diagram,  but  there  are  other  possibilities. 


In  summary,  J{b,y\u)  is  0  if  |  y-vsin(h)|  >  u.  Otherwise,  let 
5  =  -^/M^^-(y^-v^in(h^ .  Then 
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f{b,y\u)  = 


1 


4;r5(w,v) 


2vcos(/))  . 


5  + vcos(Z)) 


if  vcos{b)  >  s 


if  -5  <  vcos(/))  <  s. 


(1.14) 


0  if  vcos(Z))  <  -5 


The  conditional  density  j{y\b,u)  is  just  the  ratio  of  (1.14)  to  (1.13).  With  the  same 
exclusions  and  the  same  definition  of  5  as  in  (1.14),  it  is 


f{y\b,u)  =  < 


\l{nu)  if  vcos(Z))  >  u 
vcos(b) 


^(v cos(b)JrcCos(-v cos(b) /u)  +  ^Ju^  -v^  cos^(b)) 
.5{s  +  V  cos(b)) 


if  s<vcos{b)<u 


(1.15) 


s{v cos(b)ArcCos{-v cos{b) /u)  +  sju^  -v^  cos^(b)) 

0  if  vcos(Z))  <  -5 


if  -5  <  vcos{b)  <  s. 


Note  that  the  constant  S(u,v)  cancels  in  taking  this  ratio.  Figure  5  shows  the  conditional 
density  function  / (y  \  b,u)  for  b  =  0,  u  =  2,  and  v  >  2.  Large  values  of  Y  are  more  likely 
than  small  ones  when  b  =  0;  indeed,  the  density  / (y  |  b,u)  is  unbounded  at  the  endpoints. 
A  curiosity  is  that  fiy\b,  u)  is  a  uniform  distribution  on  the  interval  \y-u,  v+u\  when 
b  =  ±;z'/2 . 


Figure  5,  Conditional  density  of  speed  across  the  LOS,  given  detection  dead  ahead, 
for  a  target  with  speed  2,  as  long  as  ownship  speed  is  at  least  2, 
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4,5  JOINT  DENSITY  OF  B  AND  A,  GIVEN  U:  Ab^\u) 


Random  variable  X  is  the  target’s  veloeity  towards  ownship,  or  range  rate,  taken 
positive  when  the  range  is  deereasing.  Meetings  ean  only  happen  when  X  is  positive,  and 
large  values  of  X  are  potentially  alarming  to  ownship  beeause  of  the  eollision  potential, 
depending  on  the  bearing  b.  Range  rate  is  potentially  observable  through  a  Doppler  shift, 
even  when  the  range  itself  eannot  be  direetly  sensed.  The  density  that  is  derived  in  this 
seetion  is  the  first  step  in  eonditioning  on  the  observed  value  of  X. 

The  proeedure  for  determining is  similar  to  the  proeedure  for  determining 
j{b,y\u).  Since  X  is  a  function  of  both  B  and  \j/  according  to  (1.1),  we  again  have  a 
problem  in  determining  the  density  of  a  scalar  function  of  two  random  variables.  If  the 
angles  for  which  X{b,\^)=x  are  \|/i,  \|/2,  ...,  then/(h,x  |  w)  = 

where  X(\|/i)  denotes  the  derivative  of  X  with  respect  to  \|/  at  the  f'  angle.  The  required 
derivatives  are  all  X(\j/i)  =  -u  sin(\)/i).  Let  5=v  cos(h)  -  x,  so  that  u  cos(\j/,)  =  5  according 
to(l.l).  If|5|  >  u,  there  are  no  solutions  to  this  equation.  If  |5|  <  u,  there  are  two  solutions 

\|/i  and  \|/2,  for  each  of  which  |  X\y/^)  \=  -d^  .  Thus,  using  (1.12), 


f{b,x\u) 


1 

A7i:S{u,v) 


0  if  |(5'  |>  M  orx  <  0 


<  2x 


if 


d\<u  andx>0. 


(1.16) 


This  joint  density  function  is  unbounded  when  |5|  is  nearly  u,  which  corresponds 
to  extreme  situations  where  the  target  is  moving  directly  toward  or  directly  away  from 
ownship.  Figure  6  shows  a  typical  transect. 
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Figure  6:  A  graph  olf(b^\u)  versus  x  (range  rate)  for  v  =  5,  «  =  3,  and  6  =  0,  Targets 
discovered  dead  ahead  (6  =  0)  have  a  strong  tendency  to  have  a  range  rate  almost  as 
large  as  v  +  «,  with  a  somewhat  smaller  tendency  for  x  to  he  just 

above  v  -u. 


Strongly  peaked  densities,  sueh  as  those  shown  in  Figure  6,  ean  be  problematic 
for  techniques  that  evaluate  densities  only  on  a  grid.  In  some  cases,  the  cumulative  form 
of  the  density  may  be  more  useful.  It  happens  that  (1. 1 6)  can  be  integrated  analytically  on 
X  to  obtain  F{b,x\  u)  =  P{X  <x,B^db\U  =  u)/db,  a  partially  cumulative  joint 
distribution.  With  t  =  vcos(6),  this  distribution  is  given  by 


F(b,x\  u)  = 


0  ifx<0  orx<t-w, 

/(6 1  m)  if  X  >  t  +  M,  or  otherwise, 

t{7T:  -  Arc cos((x -t)l u))  - -(x-tY 
27rS(u,v) 


i{t>u 


(1.17) 


-(x-tY  +  t{Arc  cos(-t  /  w)  -  Arc  cos((x  -t)l  u))  . 


27rS{u,v) 


if  t  <  w. 


4.6  THE  PROBABILITY  OF  LEADING  THE  TARGET 

When  two  ships  meet,  assuming  that  neither  ship  subsequently  changes  either 
course  or  speed,  it  is  possible  that  the  target  will  eventually  pass  in  front  of  ownship,  at 
which  time  the  bearing  will  be  0.  In  that  case  ownship  is  said  to  be  “leading”  the  target.  If 
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the  initial  bearing  B  is  positive  (negative),  this  will  happen  if  and  only  if  Y  is  negative 
(positive);  that  is,  B  and  Y  must  have  different  signs  (all  angles  are  restrieted  to  the 
interval  [-;r,;r]  in  this  paper).  Sinee  the  event  E  of  erossing  over  one’s  bow  is  an 
important  one  for  mariners,  it  is  of  interest  to  find  its  probability.  We  have 


P(E)  =  2  j  f{b,y/\u)dbdy/ 

B>0,Y<0 


1 

2nS{u,v) 


j  (v  eos(/))  -  u  cos(^))^  dbdy/  .(1.18) 

B>0,Y<0 


The  initial  factor  of  2  in  (1.18)  is  by  symmetry — the  possibility  that  5  <  0  is  not 
accounted  for,  but  clearly  the  chances  that  the  target  will  pass  in  front  of  ownship  are  the 
same  regardless  of  the  sign  of  B. 

There  appear  to  be  only  three  cases  where  the  integral  in  (1.18)  can  be 
accomplished  analytically.  The  first  two  are  almost  obvious.  If  w  =  0  and  v  >  0,  P{E)  can 
easily  be  shown  to  be  0,  or  if  v  =  0  and  m  >  0,  P{EJ)  is  0.5.  Both  of  these  results  should 
make  intuitive  sense;  the  first  because  a  stationary  target  will  never  cross  ownship ’s  bow, 
and  the  second  because  the  direction  of  motion  for  a  stationary  ownship  should  be 
irrelevant.  The  third  case  is  when  the  two  speeds  are  equal.  One-fourth  of  such  meetings 
will  be  leading,  as  we  next  show.  Without  loss  of  generality,  we  assume  both  speeds 
are  1. 

When  u  =  v=\,  P{E)  is  MI{2tiS{\,Y)),  where  M  is  the  right-hand  integral 
in  (1.18): 

K 

M  =  ^db  ^  {cos.{b) - cos{y/))dy/ .  (1-19) 

0  Q(b) 

Here  Q{b)  is  the  set  of  angles  \\>  for  which  both  (1)  cos(/))  >  cos(\j/)  and  (2) 
sin(/))  +  sin(\j/)  <  0.  The  first  requirement  is  satisfied  when  |\|/|  >  b,  but  the  second  cannot 
be  satisfied  if  \|/  >  b  because  the  sum  of  two  positive  terms  cannot  be  negative,  so  we 
must  have  \|/  <  -b.  If  />  <  iil2,  the  second  is  satisfied  if  and  only  if  ^  e  {b-n,-b) ,  and 
Q{b)  is  therefore  that  interval.  b>  tiI2,  the  second  is  satisfied  if  and  only  if 
xj/  e  {-b,b-7i:) ,  but  being  in  that  interval  is  inconsistent  with  the  requirement  that  -b, 

ttH  -b 

so  Q{b)  is  empty.  We  therefore  have  M  =  ^  db  ^  {cos{b) -cos{x//)) di// .  It  is  now  just  a 

0  h-TT 

matter  of  evaluating  standard  integrals  to  conclude  that  M=2.  Since  >S'(1,1)  is  Ain,  we 
conclude  that  P{E)  =  0.25. 

Mariners  sometimes  refer  to  a  “lag”  situation  where  the  target’s  and  ownship ’s 
velocity  vectors  point  to  opposite  sides  of  the  LOS.  Such  situations  never  result  in  the 
target  crossing  ownship’s  bow,  so  such  events  are  not  included  in  E.  Even  if  the  target’s 
velocity  vector  points  to  the  same  side  of  the  LOS  as  ownship’s,  it  is  still  possible  that  the 
target  will  not  cross  ownship’s  bow.  Such  situations  are  sometimes  described  as 
“overlead”.  We  mention  this  only  to  emphasize  that  our  current  usage  does  not  include 
“overlead”  in  the  event  E,  which  remains  the  event  that  the  target  will  (assuming  that 
neither  ship  changes  speed  or  course)  cross  ownship’s  bow.  In  other  words,  the  opposite 
of  E  is  not  “lag”,  but  rather  “lag  or  overload”. 
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We  are  unable  to  evaluate  P{E)  analytieally,  exeept  in  the  three  eases  deseribed 
above,  but  it  is  not  diffieult  to  evaluate  numerieally.  As  expeeted,  P{E)  deereases  from 
0.5  to  0  as  ownship  beeomes  faster  and  faster  relative  to  the  target,  with  speed  equality 
being  about  “half  way  to  infinity”. 
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5. 


INFERENCES  ABOUT  R  IN  THE  CASE  OF  PASSIVE  SONAR 


Up  to  this  point,  we  have  not  been  specific  about  the  a  priori  distribution  of 
{R,U).  Other  than  the  idea  that  proximity  implies  detection,  we  have  assumed  nothing 
about  the  sensors  involved.  That  must  necessarily  change  in  this  section,  which  deals 
with  estimating  the  range  at  detection  {K).  When  detection  is  made  passively,  R  is 
notoriously  variable,  in  addition  to  being  unobservable.  Ownship  generally  becomes 
aware  of  a  target’s  presence  when  some  sound  emitted  by  the  target  is  heard  by  ownship, 
at  which  time  the  bearing  B  is  measured.  However,  the  received  signal  at  first  detection  is 
typically  near  the  threshold  of  audibility,  so  the  strength  of  the  received  signal  is  of  little 
use  in  determining  R.  Our  interest  here  is  in  making  inferences  about  R  based  on 
quantities  that  can  be  observed  at  the  time  of  detection,  particularly  B  and  its  rate  of 
change  Z.  We  therefore  must  make  a  priori  assumptions  about  R  that  are  tailored  to  the 
sensor  involved. 

The  problem  of  making  inferences  about  R  is  especially  difficult  when  detection 
is  made  by  passive  sonar,  since  the  acoustic  transmission  properties  of  the  ocean  are  hard 
to  predict.  Seemingly  one  could  reason  that  the  signal  excess,  as  predicted  by  the  sonar 
equation  (Urick,  1975),  must  be  about  0  at  the  time  of  initial  detection,  and  calculate  the 
unique  range  for  which  that  is  true.  The  problem  is  that  the  sonar  equation  itself  is  only  a 
rough  approximation  to  reality — it  can  easily  be  off  by  9  decibels  in  predicting  the 
signal-to-noise  ratio  (Wagner  et  ah,  1999,  Section  509),  which  corresponds  to  a  factor  of 
about  8  in  signal  power.  Adding  to  the  difficulty  is  that  R  depends  strongly  on  the  signal 
level  emitted  by  the  target,  which  in  turn  depends  strongly  on  the  speed  U  of  the  target, 
which  is  also  typically  unknown.  We  cannot  proceed  under  the  assumption  that  either  R 
or  U  is  a  known  quantity,  or  that  the  two  are  independent  of  each  other.  We  will  develop 
a  joint  prior  distribution  for  {R,  U)  by  first  assuming  a  distribution  for  U,  then  using  the 
sonar  equation  to  develop  a  distribution  for  R,  given  U,  and  finally  multiplying  the 
two  together. 

5.1  A  PRIOR  DISTRIBUTION  g{r,u)  FOR  PASSIVE  SONAR 

We  take  the  distribution  of  U  to  be  lognormal  with  mode  Uo;  that  is,  the  natural 
logarithm  of  (U/Uo)  is  assumed  to  be  normal  with  mean  0  and  standard  deviation  ou-  Uq 
can  be  thought  of  as  a  guess  at  the  target’s  speed,  with  quantifying  the  vagueness  of 
the  guess.  The  lognormal  distribution  is  concentrated  on  the  positive  real  line,  the  natural 
domain  of  speed,  and  has  both  scaling  (Uo)  and  shaping  (at/)  parameters. 

We  will  use  the  sonar  equation  to  develop  the  density  of  R,  given  U.  The 
signal-to-noise  ratio  depends  on  ownship  speed  v,  in  addition  to  U,  since  ownship  speed 
determines  the  level  of  noise  with  which  the  signal  must  compete.  Let  Vo  be  the  speed  at 
which  self-noise  is  equal  to  background  noise,  and  assume  that  self  noise  is  proportional 
to  the  a  power  of  ownship  speed.  Total  noise  N  is  then  proportional  to  v“+  Fo“.  The 
signal  S  is  assumed  proportional  to  l/'  for  some  parameter  P,  and  also  proportional  to  the 
inverse  square  of  R  (spherical  spreading).  The  signal-to-noise  ratio  is  then 
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(1.20) 


5  _U^(R/R,y^ 

v“+Fo“ 

where  K  is  approximately  1  and  Ro  ineludes  the  aforementioned  proportionality  eonstants. 
To  be  precise,  ln(iQ  is  assumed  to  be  normal  with  mean  0  and  standard  deviation  2gr.  If 
the  sonar  equation  has  a  normal  error  with  standard  deviation  6  decibels,  then  is  0.69. 
If  a  =  P,  can  be  thought  of  as  the  most  likely  detection  range  when  U  and  v  are  both 
equal  and  large.  We  expect  a  =  P  =  3,  since  the  power  required  to  drive  a  ship  typically 
increases  with  the  third  power  of  its  speed. 

There  is  an  important  caveat  to  the  assumption  that  K  is  approximately  1  at 
detection.  Passive  sonars  are  typically  not  omnidirectional,  with  the  location  of  the  sonar 
on  ownship  being  such  that  directions  forward  {B  =  0)  are  covered,  but  not  directions  aft 
(B  =  n).  If  the  initial  bearing  is  at  or  near  the  sonar’s  cutoff,  but  drawing  forward  (B  and  Z 
have  opposite  signs,  see  Figure  3),  then  K  might  be  substantially  greater  than  1,  and  the 
target  might  be  closer  than  would  otherwise  be  expected.  For  such  sonars,  the  Bayesian 
mathematics  developed  here  applies  only  when  the  initial  discovery  is  strictly  within  (not 
on  the  edge  of)  the  angular  coverage  region  of  the  sonar. 


Figure  7:  Illustrating  a  situation  where  the  target  is  discovered  when  entering  the 
zone  of  the  sonar’s  coverage.  The  hearing  B  is  positive  and  decreasing. 

Since  S/N=  1  at  first  detection,  (1.20)  implicitly  defines  the  density  of  R,  given 

U  =  u.  Specifically,  ln(i?/i?o)  is  normal  with  mean  0.5  ln( - )  and  standard  deviation 

v“  +  Fo“ 

Or.  Parameters  that  need  to  be  known  for  determining  the  joint  distribution  of  R  and  U 
are  thus  Uq,  Vo,  Ro,  ol,  P,  cju,  and  or.  Given  those  parameters,  the  joint  density  function  is 

c::p(  Jln(r/.S„)-0.51n(«'’/(v'-+K„'-))f,^ 

g(r,u)  = - ^ ^ - .  (1.21) 

InrucT^jCTg^ 

This  density  function  is  illustrated  on  sheet  “Prior”  of  Meetings.xls,  with 
opportunities  to  explore  the  impact  of  different  parameter  choices. 
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There  are  other  unimodal  distributions  with  the  appropriate  scaling  and  shape 
parameters,  but  (1.21)  uses  standard  functions  and  has  some  convenient  analytic 
properties.  One  advantage  of  the  lognormal  assumption  is  that  the  function  R^u)  defined 
in  Section  1.2  can  be  expressed  analytically  (see  Johnson  &  Kotz,  1970),  as  can  certain 
other  moments.  Some  useful  formulas  are  recorded  below: 


E{U)  =  U^qxy){.5ct^j') 


^(^)  =  ^0 .  exp(0.5rT/ +  0. 125(y0r7^)^) 

V  V  r  f. 


^^^exp(2<T/  +0.5(/?<t„)") 

V  +  r  Q 


E{RU)  =  R,U  +0.5{{l  +  0.5/3)a^f) 

y  V  +  Kq 


Rgin)  =  \rg(r,u)dr 


R. 


2  2. 


exp(0.5[rT/-(^ - ^)^]). 


cr„ 


(1.22) 

(1.23) 

(1.24) 

(1.25) 

(1.26) 


It  may  be  tempting  to  use  the  experience  of  subject  matter  experts  (SMEs)  to 
estimate  some  of  the  parameters  required  for  this  sonar  model,  but  a  word  of  caution  is 
appropriate.  It  is  likely  that  SMEs  will  have  more  knowledge  of  targets  that  are  met  than 
of  targets  in  general,  but  the  formulas  in  this  section  are  for  targets  in  general.  It  would  be 
a  mistake,  for  example,  to  use  (1.22)  with  E(f/)  =  5  kt  because  of  an  SME’s  statement 
that  “targets  I  have  met  have  an  average  speed  of  5  kt”.  Such  SME  knowledge  applies  to 
distributions  upon  meeting,  the  subject  of  the  next  section,  rather  than  the  a  priori 
distributions  that  are  the  subject  of  this  section. 

5.2  DISTRIBUTIONS  UPON  MEETING 

The  joint  density y(^,M)  upon  meeting  is  given  by  (1.8),  with  (1.21)  substituted  for 
g{r,u).  We  know  of  no  analytic  expression  for  the  constant  A  from  (1.8),  so  are  forced  to 
resort  to  numerical  integration,  but  otherwise  this  density  is  an  analytic  function.  Various 
moments  of  R  and  U  can  be  derived  by  numerical  integration,  and  possibly  used  for 
parameter  estimation  through  SME  judgments. 

The  joint  density  Xem, 6)  is  the  product  of  (1.8)  and  (1.13).  The  equivalent  speed 
S{u,v)  cancels  when  making  that  product,  so  J[r,u,b)  is  again  an  analytic  function  except 
for  A.  If  the  bearing  is  observed  to  be  b  upon  meeting,  then  this  joint  density  easily 
enables  the  computation  of  the  conditional  density  y(r,M| 6)  and  its  marginals  1 6)  and 
J{u\b).  In  other  words,  everything  that  is  to  be  known  about  the  range  and  speed  of  the 
intercepted  target,  given  only  the  bearing  to  the  target,  is  at  hand. 
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The  joint  dQmiiy  j{r,u,b,y)  is  the  produet  of  (1.8)  and  (1.14).  Let  the  bearing  rate 
be  Z,  so  that  Z=YIR.  The  joint  density  y(r,M,h,z)  is  then  rj{r,u,b,y  =  rz),  where  this 
awkward  notation  means  that  rz  is  to  be  substituted  for  y.  We  ean  now  eondition  on  (b,z) 
to  diseover  the  effeets  of  observations  of  bearing  and  bearing  rate  on  the  distribution  of 
iR,U).  This  is  the  “instant  TMA”  referred  to  in  the  introduetion. 

The  joint  density  J{r,u,b,x)  is  the  produet  of  (1.8)  and  (1.16),  enabling  the 
eomputation  of  f{r,u\b,x),  the  density  of  range  and  speed  when  both  bearing  and  range 
rate  are  observed. 

The  distributions  deseribed  above  are  all  implemented  and  exereised  in  the 
workbook  Meetings.xls.  See  Seetion  7. 
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6.  VISUAL  BASIC  FOR  APPLICATIONS  (VBA) 
IMPLEMENTATION  OF  VARIOUS  DISTRIBUTION  FUNCTIONS 


The  VBA  module  ineluded  with  RandMeet.xls  eontains  eode  for  all  of  the 
funetions  introdueed  above,  plus  some  additional  ones  useful  for  simulation.  In  that 
module,  densities  defined  upon  meeting  all  start  with  the  letter  f,  followed  by  the 
upperease  names  of  the  random  variables  whose  joint  density  is  under  eonsideration, 
possibly  followed  by  the  letter  g  (for  “given”),  followed  by  the  upperease  names  of  the 
random  variables  whose  values  are  given.  Thus,  fXgY(x,y)  is  the  name  of  the  funetion 
that  would  be  ealled/yy(x[y),  were  subseripts  available  in  VBA  eode.  Some  densities  omit 
normalization  by  hard-to-eompute  eonstants  sueh  as  A,  in  whieh  ease  the  name  of  the 
funetion  ends  in  the  letter  k  as  a  warning. 

6.1  SECTION  2  FUNCTIONS 

S{u,v)  (formula  1.4)  is  EquivSpd(M,v). 

6.2  SECTION  4  FUNCTIONS 

f{b\u)  (formula  1.13)  is  fBgUk(Z),M,v).  This  funetion  needs  to  be  divided  by 
2S(m,v)  to  be  a  true  density.  The  alternative  funetion  BearDen(Z),M,v)  is 
Koopman’s  original  speeifioation,  and  is  a  true  density. 

f{\\>\b,u)  (the  ratio  of  formula  1.12  to  1.13)  is  fPgBU(M,v,/),  p),  a  true  density. 

f{y\b,u)  (formula  1.15)  is  lY gB\J{u,v,b,y),  a  true  density. 

A^a|w)  (formula  1.16)  is  fBXgUk(M,v,6,x).  Divide  by  27t:S(m,v)  to  get  a 
true  density. 

F(b,x\u)  (formula  1.17)  is  fBXgUo(u,v,b,x).  The  terminal  “e”  means 
“eumulative”.  Divide  by  27rS(M,v)  to  get  the  true  distribution. 

g{r,u)  (formula  1.21)  is  prior(r,  RO,  u,  UO,  v,  VO,  SigU,  SigR,  j!70w).  Parameters  a 
and  P  are  both  assumed  to  be  equal  to  input  pow.  The  argument  SigU  is  a^,  SigR 
is  <jR,  and  other  symbols  are  as  used  earlier. 

E(U)  (formula  1.22)  is  EU(R0,  UO,  v,  VO,  sigU,  SigR,  pow),  notation  as  in  g{r,u). 
E(R)  (formula  1.23)  is  ER(R0,  UO,  v,  VO,  sigU,  SigR,  pow),  notation  as  in  g(r,u). 
E(R^)  (formula  1.24)  is  ER2(R0,  UO,  v,  VO,  sigU,  SigR,  pow),  notation  as  in  g(r,u). 
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E(RU)  (formula  1.25)  is  ERU(i?0,  UO,  v,  VO,  sigU,  SigR,  pow),  notation  as 
mg{r,u). 


Rg{u)  (formula  1.26)  is  Rg(i?0,  u,  UO,  v,  VO,  sigU,  SigR,  pow),  notation  as 
mg{r,u). 

A  (formula  1.6  ,  using  Rg{u)  and  S(m„v)  from  above)  can  be  obtained  using  the 
function  called  mtegral(v,  RO,  SigR,  UO,  sigU,  pow,  VO),  with  notation  as  in 
g{r,u).  This  is  a  numerical  integral,  so  comparatively  expensive  to  evaluate. 

E(R|meetmg)  (not  numbered)  can  be  obtained  from  ERA(v,  RO,  SigR,  UO,  sigEi, 
pow,  VO),  notation  as  in  g{r,u).  This  numerical  integral  is  the  expected  range  of  a 
met  target  when  ownship  moves  at  speed  v,  and  might  be  of  use  in 
parameter  estimation. 

j{r,u,b)  (not  numbered)  is  fREiBk(r,  RO,  u,  UO,  v,  VO,  sigU,  SigR,  pow,  b).  This 
needs  to  be  multiplied  by  47rac/Ji;A  to  be  a  true  joint  density. 

j{r,u,b,z)  (not  numbered)  is  fREiBZk(r,  RO,  u,  UO,  v,  VO,  sigU,  SigR,  pow,  b,z). 

6.3  SAMPLING  AND  OTHER  FUNCTIONS 

Various  other  functions  are  also  provided  in  Meetings. xls,  generally  introduced  by 
comments  in  the  VBA  code.  Among  these  are  some  sampling  functions  that  may  be  of 
use  to  those  who  have  no  ambition  to  do  any  Bayesian  analysis,  but  would  like  to 
generate  random  meeting  parameters  in  a  Monte  Carlo  simulation.  All  use  the  rejection 
method  (Ross,  2000).  The  functions  are: 

Eor  B  given  U,  use  Bsample(M,v)  or  possibly  Bsample2(M,v).  The  latter  should  be 
somewhat  faster,  but  the  code  is  more  complicated. 

Eor  C  given  U  andB,  use  Csample(M,v,h).  If  {u,v,b)  is  an  impossible  combination, 
this  function  will  return  a  very  large  positive  number. 

Eor  Y  given  U  and  B,  use  Ysample(M,v,h).  If  (u,v,b)  is  impossible,  this  function 
will  return  a  very  large  positive  number. 

Eor  both  B  and  C  given  U,  use  subroutine  BCSample(M,v,6,c). 

These  functions  are  not  exhaustive,  but  others  can  be  based  on  densities  derived 
above.  Eor  example,  (1.7)  could  be  used  to  give  a  joint  sample  of  {R,B,U,C). 
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7. 


THE  WORKBOOK  MEETINGS.XLS 


This  workbook  includes  several  sheets  named  after  density  funetions  deseribed  in 
Seetion  6.  Eaeh  of  these  sheets  is  a  graphieal  or  eomputational  exploration  of  the  funetion 
with  its  name.  Sheet  “leading”  permits  the  approximate  eomputation  of  P{E)  as  deseribed 
in  Seetion  4.6.  These  are  all  “Seetion  4”  sheets;  that  is,  no  referenee  is  made  to  a  prior 
distribution  for  {R,U). 

The  five  “range”  sheets  whose  names  start  with  R  are  all  eontrolled  by  the  sheet 
named  “RUBZ”,  whieh  ineludes  a  eommand  button  that  updates  all  five.  The  range 
sheets  empower  the  exploration  of  the  extent  to  whieh  inferenees  about  R  ean  be  made, 
based  on  observations  of  B  and  possibly  X  or  Z.  Speoifieally, 

sheet  RUprior  shows  the  prior  density  of  {R,U)  before  meeting; 

sheet  RU  shows  the  density  for  met  targets; 

sheet  RUB  shows  the  density  eonditioned  on  knowing  5; 

sheet  RUBX  shows  the  density  eonditioned  on  knowing  B  and  X;  and 

sheet  RUBZ  shows  the  density  eonditioned  on  knowing  B  and  Z. 

The  joint  distribution  of  {R,U)  is,  in  all  oases,  displayed  digitally  after  multiplying 
by  a  large  oonstant  (100)  to  improve  readability.  Numbers  larger  than  0.1  are  oolored  to 
get  a  graphieal  idea  of  the  region  in  whieh  density  is  “signifioant”.  In  all  oases,  the 
marginal  distribution  of  R  is  oomputed  by  summing,  then  used  to  approximate  E(R)  under 
the  speoified  oondition,  and  then  oopied  to  the  oontrolling  sheet  RUBZ  to  permit  easy 
oomparisons.  In  some  oases,  the  exaot  value  of  E{R)  is  also  known  and  used 
for  oomparison. 

The  range  pages  have  many  inputs,  so  many  experiments  are  possible.  One 
experiment  begins  with  the  “Standard”  soenario  (Tools  menu),  exploring  the 
oonsequenoes  of  modifying  oonditions  on  B  and  Z.  The  a  priori  mean  deteotion  range 
upon  meeting  in  this  soenario  is  1.52  nm.  Table  1  shows  the  mean  deteotion  range  with 
various  observations. 


Bearing  (B] 

1  in  radians 

Bearing  rate  (Z)  in 
radians  per  hour 

0 

0.5 

1 

1.5 

Null 

1.40 

1.41 

1.48 

1.67 

0 

2.03 

1.72 

2.75 

3.35 

1 

2.09 

1.88 

2.36 

2.80 

2 

1.50 

1.74 

2.13 

2.34 

3 

4.27 

1.40 

1.76 

1.84 

4 

0.77 

1.12 

1.45 

1.49 

Table  1:  The  mean  range  in  nm  for  several  assumptions  about  bearing  and  bearing 
rate,  “Null”  meaning  that  no  observation  of  Z  is  made. 
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Table  1  displays  an  anomaly  when  Z  =  3  and  5  =  0,  where  the  expected  range  is 
clearly  not  related  to  the  values  in  neighboring  cells.  The  anomaly  is  easily  explained,  but 
nonetheless  is  something  to  beware  of  in  computations  such  as  these  involving 
unbounded  densities  on  a  grid,  which  in  this  case  includes  values  of  R  and  U  that  are 
multiples  of  0.1.  When  Z  =  3  and  5  =  0,  it  happens  that  the  point  (R,U)  =  (4.3,  12.9)  is 
nearly  on  a  singularity  of  the  joint  density,  to  the  extent  that  it  ends  up  being  the  only 
possibility  with  any  significant  probability.  This  one  point  takes  over  the  numerical 
computation  of  the  mean  value  of  R,  forcing  it  to  be  about  4.3.  An  operational  version  of 
these  computations  would  recognize  and  protect  against  such  anomalies,  but  our  simple 
grid-based  computation  does  not. 

Except  for  the  anomaly.  Table  1  shows  a  decrease  in  expected  range  as  the 
bearing  rate  increases.  This  should  be  intuitive  because  Z  is  inversely  proportional  to  R. 
However,  the  decrease  is  not  as  fast  as  one  might  expect,  and  is  certainly  not  inversely 
proportional  to  Z.  The  reason  for  this  mushiness  is  that  there  are  other  explanations  for  a 
high  bearing  rate:  the  target  might  be  going  fast,  or  the  target’s  relative  course  \|/  might  be 
nearly  n/2.  Since  neither  U  nor  \|/  has  been  observed,  it  makes  intuitive  sense  to  lend 
some  credence  to  both  of  those  possibilities,  in  addition  to  the  possibility  that  R  might 
be  small. 

Table  1  (the  Null  row)  also  shows  that  the  expected  range  increases  with  the 
observed  bearing.  Intuitive  or  not,  that  is  the  truth  of  the  matter.  Although  not  shown  in 
Table  1,  each  of  the  quantities  displayed  there  is  surrounded  by  a  widespread  distribution. 
Figure  1  displays  three  distributions  associated  with  the  Standard  scenario,  all  of  which 
exhibit  significant  variance. 

The  results  in  this  section  are  typical,  and  demonstrate  the  difficulty  that  ownship 
must  inevitably  have  when  meeting  a  new  target,  prior  to  making  any  maneuver.  Instant 
TMA  will  not  normally  be  definitive,  so  considerable  risk  will  remain  if  decisions  must 
be  made  quickly.  To  accurately  predict  the  target’s  future  movement,  it  will  be  necessary 
to  maneuver  and  accept  the  implied  time  delay. 

Probability  density  of  range 


-  -  —  Detect 

"Bearing&x” 


0  1  2  3  4  5  6 


Figure  8:  Marginal  densities  of  range  under  three  variations  of  the  Standard 
scenario.  The  “Detect”  curve  assumes  only  detection,  and  the  “Bearing&x”  curve  is 
conditioned  on  knowing  the  hearing  and  range  rate.  All  three  exhibit 

significant  variance. 
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